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ABSTRACT 

Aims. Powerful gravitational telescopes in the form of massive galaxy clusters can be used to enhance the light collecting power over 
a limited field of view by about an order of magnitude in flux. This effect is exploited here to increase the depth of a survey for lensed 
supernovae at near-IR wavelengths. 

Methods. We present a pilot supernova search programme conducted with the ISAAC camera at VUT. Lensed galaxies behind the 
massive clusters A1689, A1835, and AC114 were observed for a total of 20 hours divided into 2, 3, and 4 epochs respectively, 
separated by approximately one month to a limiting magnitude y<24 (Vega). Image subtractions including another 20 hours worth of 
archival ISAAC/VLT data were used to search for transients with lightcurve properties consistent with redshifted supernovae, both in 
the new and reference data. 

Results. The feasibility of finding lensed supernovae in our survey was investigated using synthetic lightcurves of supernovae and 
several models of the volumetric Type la and core-collapse supernova rates as a function of redshift. We also estimate the number 
of supernova discoveries expected from the inferred star-formation rate in the observed galaxies. The methods consistently predict 
a Poisson mean value for the expected number of supernovae in the survey of between Nsn=0.8 and 1.6 for all supernova types, 
evenly distributed between core collapse and Type la supernovae. One transient object was found behind A 1689, 0.5" from a galaxy 
with photometric redshift Zgai = 0.6 ± 0.15. The lightcurve and colors of the transient are consistent with being a reddened Type IIP 
supernova at zsn = 0.59. The lensing model predicts 1.4 magnitudes of magnification at the location of the transient, without which 
this object would not have been detected in the near-IR ground-based search described in this paper (unlensed magnitude J ~ 25). 
We perform a feasibility study of the potential for lensed supernovae discoveries with larger and deeper surveys and conclude that 
the use of gravitational telescopes is a very exciting path for new discoveries. For example, a monthly rolling supernova search of a 
single very massive cluster with the HAWK-I camera at VLT would yield >10 lensed supernova lightcurves per year, where Type la 
supernovae would constitute about half of the expected sample. 

Key words, cosmology: gravitational lensing supernovae: general 



1. Introduction 

Acting as powerful gravitational telescopes, massive galaxy 
clusters o ffer unique opportu nities to observe extremely distant 
galaxies dKneib et al.L l2004l) . as well as distant supernovae 
(SNe) , t oo faint to be otherwi se detected ( Kolatt & Bartelrnani] , 



19981 iKovner & Paczvns kil 119881: IS ullivan et all I200C; 



Gal- Yam et al.L 12002: Gun narsson & Goob ar. 2003^ Lensing 
magnifications of up to a factor ~40 have been inferre d for 
many multiple images of galaxies (ISeitz et all Il998l) and 
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typical magnification factors of 5 to 10 are common within the 
central few arcminutes of the most massive clusters of galaxies. 
Exploiting this remarkable boost in flux is an interesting avenue 
for probing the rate of exploding stars at redshifts beyond 
the detection capabilities of currently available telescopes. 
Successful programs detecting intermedi ate to hi gh reds hift SN 
at optic al wavelengt hs include SDSS-II (iFriem an et al.. l2008l). 
SNLS (lAstier et al.L l2006h . and ESSENCE dMiknaitis et al.L 
|2007|) . which al l target zS^ super novae, mainly Type la. For 
higher redshifts, iRiess et aTl ( l2007h demonstrated the power of 
space measurements by reporting the discovery and analysis of 
23 Type la SNe with z > 1, although not without significant 
effort. The project made use of about 750 HST orbits to detect 
SNe, obtain multi-color lightcurves and grism spectroscopy 
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with ACS and NICMOS. Dawson et flL(in prep.) improved 
the yield of high-z Type la SN detections with ACS/HST by 
targeting SN in massive z > 1 clusters. A common feature of 
both HST projects was that the search was done in the F850LP 
filter, i.e., z-band. iPoznanski et alj (|2007|) used the extremely 
large FoV of the Suprime-Cam at Subaru 8.2m to measure SNIa 
rates up z = 1.6. They reported the discovery of 13 SNla beyond 
z = 1 with repeated imaging of the Subaru Deep Field at optical 
wavelengths, including z-band. 

We explore a complementary technique for SN detection and 
photometric follow-up that involves a near-lR, "rolling", SN sur- 
vey behind intermediate-redshift massive clusters. By exploit- 
ing the significant lensing magnification, the redshift discovery 
limit is enhanced. There are, however, two main limitations to 
this approach. First, the large lensing magnification is limited 
to small solid angles around the cluster core, of typically a few 
arc-minutes. Second, conservation of flux implies that the sur- 
vey area behind the cluster in the source plane is shrunk due to 
lensing. 

The choice of near-lR filter ensures that the survey has a 
potential for SNla discoveries to unprecedented high redshifts, 
z ~ 2.5, and still samples the rest-frame optical part of the spec- 
trum, potentially allowing a significant increase in the lever arm 
of the Type la SN Hubble diagram. Furthermore, because of 
the strong lensing effect, multiple images of high-z SNe with 
time separations of between weeks and a few years could be 
observed. These rare events could provide strong const raints on 
the H ubble constant, using the time-delay technique dRefsdal 
1 19641) and possibly be used as tests of d ark matter and dark 
energy in an unexplored r edshift range (iGoobaretali l2002t 
iMortsell & Sunessonl l2006h . A feasibility study of the poten- 
tial for improving the mass models of clusters of galaxies us- 
ing lensed SNe will be presented in Riehm et ai, in preparation 
(Paper 111). 

In this paper, we describe a pilot program using the ISAAC 
near-lR imaging camera at ESO's Very Large Telescope (VLT), 
to detect gravitationally lensed SNe behind very massive clusters 
of galaxies. The potential for a scaled-up version of this project 
with the new HAWK-1 near-lR instrument at VLT is also stud- 
ied. A description of the observing strategy, the data-set and data 
reduction, as well as a full presentation of the photometry of the 
transient object dis cussed in Sect. |6.2[ ar e presented in an ac- 
companying paper dStanishev et al.l (l2009l) . Paper 1). 

Throughout this paper, we adopt the concordance model cos- 
mology, Q.M = 0.3, Qa - 0.7, li = 0.7, w = -1. Magnitudes are 
given in the Vega system. 



2. Supernova subtypes 

SN explosions are broadly divided into two classes, core- 
collapse supernovae (SN CC), marking the end of very massive 
stars (>8Mo), and the so-called Type la supernovae (SNla), be- 
lieved to be either the result of merging white dwarfs or accreting 
white dwarfs in close binaries, where thermonuclear explosions 
are triggered when the system is close to the Chandrasekhar 
mass, 1.38M0. 

Several subtypes of explosions belong to the core-collapse 
class, including Type 11 SN as well as Type Ib/c and hypernovae 
(HN), which are also interesting because of their association 
with GRBs. Type 11 SNe are furthermore subdivided into IIP, 
IIL, and lln based on lightcurve and spectroscopic properties. 
For reviews of SN clas sification and their general properties, see 
lFilippenkd(ll997h and lLeibundgu3(l2008h . 



Table 1. Supernova properties. 



SN type 


My (mag) 


a-My (mag) 


fee 


la 


-19.23 


0.30 




np 


-16.90 


1.12 


0.50 


nL 


-17.46 


0.38 


0.2025 


HLbright 


-19.17 


0.51 


0.0675 


Iln 


-19.05 


0.50 


0.05 


Ib/c 


-17.51 


0.74 


0.15 


HN 


-19.20 


0.30 


0.03 



In Table [T] we summarize some of the main properties of the 
SNe being considered in this analysis, which are: the peak V- 
band brightness. My', the one-standard-deviation range around 
the peak intrinsic luminosity, cr^^, (a Gaussian distribution has 
been assumed for all types, except for Type IIL supernovae, 
for which a bi-Gaussian distribution is used, with the two peak 
values being labeled IIL and llLbiight ); and the fraction of the 
core-collapse SN subtypes, fee, inferred from measurements of 
the local universe. We adopted the values of fee, lightcurve 
and spectral properties con ipiled by Peter NugenQ , which in 
turn are based on work by iRichardson et al.l (|2002|) . We note, 
however, that the uncertainty in the relative frac tions within 
the co re-collapse types is quite large. For instance, ISmartt et al.l 
(I2009h found a much h igher fraction of Type Ib/c (29%) than 
IRichardson et al.l (|2002|) . while their estimate of the number of 
type IIL is about a factor 10 lower than what has been assumed 
here. Clearly, significantly larger data-sets are needed to deter- 
mine the CC rates accurately, both at low and high redshifts. 

Since the peak brightness, lighcurve shape and spectral en- 
ergy density vary significantly between SN types, we treat 
each subtype separately when computing the expected rates. 
Synthetic lightcurves of SNe at a luminosity distance diiz) 
are calculat ed for the observ er NIR filters using cross-filter K- 
corrections (iKim et all 1 19961) 



niyiz, t) = Mv(t) + D(z) + Kvy(z, t) + Am(z), 



(1) 



where Y represents an arbitrary observer filter and the distance 
modulus is defined as 



£)(z) = 25 + 51ogi„ 



ddz) 
IMpc 



(2) 



The luminosity distance, d^, in a flat (Q.^ - 0) Friedmann- 
Lemaitre-Robertson- Walker model of the universe is given by: 



dL(z)^c{l+z) 



r dz' 

Jo Hiz'Y 



(3) 



where c denotes the speed of light in vacuum and the Hubble 
parameter evolves with redshift as 



H{z) = 100 ■ h VQm(1 +z)UQa , 



(4) 



in units of km s^'Mpc"'. 

We also included in Eq. (Tj the perturbation Am(z) = 
Amext(z) + Amiens(z) in the observed magnitude from lensing 
magnification, Amiens, and/or extinction by dust along the line 

of sight, AWext- 

Figure[T]shows examples of synthetic lightcurves for a repre- 
sentative set of SN types at redshifts z = 1.5 and z = 2.0 through 
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Fig. 1. Synthetic lightcurves in S'Z-band of Type la. Type Ib/c, 
and Type IIP supernovae at z = 1 .5 and z = 2 at their mean peak 
brightness. Also shown is the increased sensitivity for lensing 
magnification by 1, 2, 3 and 4 mag. The unlensed limit, ~ 24, 
corresponds to the search depth for most of the data in Table [3] 
Type IIP supernovae are too faint to be detected with the ISAAC 
survey at z = 2, even with Amiens - -4. 



the ISAAC 5Z-filter (A^fj = 1.06 yum; FWHM= 0.13 yum), a 
broader and redder version of the more common Y filter Due 
to the large lensing magnifications from the foreground cluster 
(Sect.lJTJ, significantly higher redshifts can be probed, not only 
for the intrinsically fainter core-collapse supernovae, but also for 
Type la supernovae beyond z ~ 1.5. 



3. Supernova rates 

We now consider two different routes for computing the ex- 
pected number of SNe in a survey. In Sect. 13.11 the volumetric 
approach is followed, i.e., the predictions are derived from the 
volume probed in the field of view of the survey and assumptions 
about the SN rate per co-moving volume for the various types 
of SNe as a function of redshift. In Sect. 13.21 we also consider 
the rates derived from the rest-frame UV luminosities of the re- 
solved galaxies behind the clusters based on the assumption that 
they trace the star-formation rate in each individual galaxy. 

3.1. Volumetric rate estimate 

The expected number of SNe of a certain subclass, dNj, in a red- 
shift interval, dz, depends on the monitoring time for that specific 
SN type, Tj, the solid angle of the survey, co, and the volumetric 
SN rate, (with units Mpc^-'yr"'), given by 



dNj = Tjiz) ■ 



(l+z) 



Furthermore, it is a function of cosmological parameters, since 
it includes the comoving volume element 



dVr 



cdliz) 



H{z){\+zY 



tOJ dz. 



(6) 



Next, we explore the current estimates of the volumetric rates 
of core-collapse and Type la supernovae. 

3.1.1. Core-collapse SNe 

Large scale SN programs such as SDSS, SNLS, ESSENCE, or 
GOODS/PANS and even the planned survey at the LSST are 
rather inefficient at detecting core-collapse SNe at z > 1. As 
an example, the two highest-redshift identified CC SNe in the 
five-year SNLS survey are at z = 0.617 (a probable Ib/c) and 
z - 0.605 (Ib/c confirmedfl We note, however, that SNLS 
specifically targeted Type la supernovae, and thus not optimized 
for finding CC SNe. 

The magnification provided by foreground clusters could en- 
able the exploration of this population for the first time. Since 
the progenitors of CC SNe are massive short-lived stars, the CC 
SN rate, r^*^, reflects the ongoing star-formation rate (SFR, units 
M0yr 'Mpc^^*). Thus, we can use the SNR to obtain indepen- 
dent bounds on the cosmic SFR since 



'•v^(z) = ^8° X SFR(z), 



(7) 



where k^^ — 0.007Mq' is estimated using a Salpeter IMF 
(ISalpeteit[T955h and a progenitor mass r ange of between 8 and 
50 solar masses, as in lDahlen et al.l(l2004l) . Although straightfor- 
ward in principle, large uncertainties plague the procedure out- 
lined in Eq. (|7|. The e stimates of SFR(z) fro m various data sets 
show a large span ([ Charv & Elbaz, 2001 ; Giaval isco et aUl2004t 
iHopkins & Beacoin , ,2006; .Mannucci et al.., ,20071) . thus leading 
to a very uncertain range of predictions for the SN rates, as 
shown in Fig.|2l but also allowing for the possibility to constrain 
the SF R(z) by measuring r^'-. The estimate by Mannucci et al] 
(|2007|) (M07) shown in Fig|2] incorporates a strong efficiency 
cut due to dust obscuration. The underlying assumption is that 
star formation correlates with dust density. Thus, as the star for- 
mation increases with redshift, the fraction of obscured SNe in- 
creases. M07 parametrized the fraction of observable CC SNe at 
optical wavelengths as being / = 0.95 - 0.28 x z for z < 2. In a 
similar way to the approach followed bv iLien & Fields! (120091) . 
the CC rate in M07 is extended smoothly up to z ~ 4 in a man- 
ner compatible with the upper limits on the fraction of radiation 
escaping high- redshift galaxies, / ~ 0.02 for the redshift inter- 
val 3 < z < 9 jGnedinet al.Ll2008h . While this is a rather con- 
servative approach, given that these estimates were derived for 
considerably shorter rest-frame wavelengths than those relevant 
to our work, it has a negligible impact on our results. 

The observation al results for the CC rates at high-z 
(iDahlen et al.Ll2004l) show an increase in the range z ~ 0.3—0.7, 
which is consistent with independent estimates of the SFR. 
Extending these measurements to z > 1 is clearly important for 
checking the underlying SFR models. Observations at near-IR, 
i.e. in the rest-frame optical, should also provide a direct probe 
of the star formation that could be missed by UV surveys due to 
extinction by dust along the line of sight. Since the VLT/ISAAC 
survey had very limited sensitivity beyond z > 2, the smoothly 
extrapolated model of M07 is used mainly for our estimate of the 
feasibility of discovering lensed SNe in future surveys in Sect.|2] 



■dVc 



(5) 
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Fig. 2. Predictions for the CC SN rate, r'^^iz) derived from var- 
ious estimates of the star formation rate, SFR(z). The shaded 
region (M07) is an ext r apolat ion based on the dust corrected 
rate in iMannucci et al ] (l2007h . The dashed line (D04) gives 
the best fit to the GOODS CC discoveries, as shown in 
iDahlen et aL (2004). T he dot t ed Un e (KN08) shows the predic- 
tion in lKobayashi & Nomotol (l2008l) . 



3.1.2. Type la SNe 

While SNIa have been used extensively for deriving cosmolog- 
ical parameters, it is unsatisfactory that the progenitor scenario 
preceding the SN explosion is still mostly unknown. Existing 
models predict that different scenarios, such as the single de- 
generate and the double degenerate models, should have differ- 
ent delay-time distributions 0(0, describing the time between the 
formatio n of the progen itor star and the ex plosion of t he SN . The 
results of lMannucci et al. (2005) and Sulli van et al.l (|2006|) sug- 
gest that the specific SNR (SNR per unit mass) is significantly 
higher in young sta r forming galaxies th an in older galaxies . 
iPritchet et aTl dioOSl) used the results from lSulUvan et alj ( |2006|) 
to show that the delay-time distribution, 0(f), is consistent with a 
power-law function and that the specific SNR is at least a factor 
of 10 higher in active star-forming galaxies compared to passive 

galaxies. 

Using a different method, IStrolger et alj ( |2004|) compare the 
SFR(t) and the SNR, r(;'(f), derived in the GOODS fields to de- 
rive the delay-time distribution with the relation 



Redshift 



Jo 



SFR(f')0(f-f')df', 



(8) 



where y is the number of SNe per unit stellar mass formed. 
Assuming that 0(f) has a Gaussian shape, they found a preferred 
delay time of r ~3-4 Gyr, which is significantly longer than that 
found when using the specific SNR described above. 

We note that the main drive r of the relatively long delay-time 
found in the method used by IStrolger et al.l ( |2004|) is the low 
number of Type la SNe found at hi gh re dshift z > 1.4. Using 
the extended GOODS survey, Dahlen et al. (2008) also found 
fewer Type la SNe at z > 1 .4 than what is expected if the delay- 
time is short and SNR follows the S FR. Results for t he Typ e la 
rate at z >1.4 were also p resented in lPoznanski et al.l (|2007|) and 
iKuznetsova et alj ( l2008h . While both these sets of results show 
a rate that is consistent with being constant at z > 1 .4, the large 
statistical errors can not e xclude a rate t hat de clines sharply as 
suggested by the results in lDahlen et al.l (l2008h . 



Fig. 3. Extrapolations of available SNIa rate predictions, rj?(z). 
The shaded regi on (M07) is an extrapo lation based on the dust 
corrected rate in IMannucci et al. I (120071) . The das hed line (DOS) 
gives the best fit to the HST supernova survey (iDahlen et al. . 
2008) corresponding to r = 3.4 Gyr The dot-dash line (TAUl) 
shows the corresponding rate for t = 1.0 Gyr. In both cases, 
a Gaussian distribution was assumed for t. The thic k solid line 
(N06) corresponds to the SNLS published rates in iNeill et all 
( 20061) . The thin sohd line (SB05) s tems from the so-called 
"A-hB model" in IScannapieco & Bilds ten (20 (31) and the dot- 
ted lin e (KN08) shows the prediction in iKobayashi & Nomotol 
(120081) . 



It is therefore particularly important to search for SNe at 
these redshifts, where the predicted rate is most sensitive to the 
delay time. If t is large (>3-4 Gyr), there should be a steep de- 
cline in the SNIa rate at z >1.5, while if r is short, the SNR 
should follow the SFR and remain fairly constant to z>3. 

Furthermore, there are also theoretical predictions that the 
SNIa rate could be significantly suppressed (or even inhibited) 
at high redshifts (z > 2 in spirals and z > 2.5 in ellipti - 
cals) due to metallicity effects, e.g., IKobayashi et al.l (Il998h . 
Kobayashi & Nomoto (2008) revised their analysis and found 
an expected increase of the SNIa rate in elliptical hosts above 
z = 2.5. These results also show that deriving the Type la rate at 
high redshift is of great interest. 

To calculate the number of detectable Type la SNe, we 
use a sample of SNIa rate model predictions and best-fit solu- 
tions to the available data, extrapolated to very high redshifts. 
These models are shown in Fig. [3] As for CC SNe, we used the 
smoothly extrapolated M07 model as a benchmark for the feasi- 
bility studies in Sect. [T] 

3.2. SN rates derived from the SFR of observed galaxies 

Since the rate of SNe is expected to follow the star-formation 
rate, we also consider the numbers that can be derived for the 
SFR in the galaxies detected along the field of view. In Sect. jSj 
we describe how galaxy catalogs were generated for the resolved 
objects along the Une of sight to massive clusters. 

We used the rest-frame UV luminosity as a tracer of the SFR 
in the observed galaxies, redshifted to the optical bands. Since 
the UV luminosity is dominated by the most short-lived stars, it 
is closely related to star formation. 
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We use L2800, the flux at rest-frame A^jf - 2800A, to esti- 
mate the SFR. We first used the photometric (or spectroscopic) 
redshift (see Sect. |5]l to derive which two observed filters strad- 
dle the rest-frame L2800 and interpolate between those using the 
best-fit spectral template to derive the apparent magnitude cor- 
responding to the rest-frame L2800 flux. The absolute L2800 mag- 
nitude is thereafter derived after correcting for distance modulus 
and K-corrections. Next, the lensing magnification is taken into 
account, as described in Sect. 14.11 

Finally, we use the relation between L2800 and SFR from 
dPahlen et al.L |2007|) to relate the flux to star formation. 



SFR(Moyr-') = 



^28Qo(erg-s 'Hz ) 
7.0 ■ 1027 



(9) 



The expected number of core collapse SNe are calculated 
using Eq. (|7]), whereas the SNIa rate is estimated using Eq. ([8j 
for a Gaussian distribution of (^(f) with (t, cr) — (3.4,0.68) in 
units of Gyr. 



4. Clusters as gravitational telescopes 

We have investigated the use of some of the most massive clus- 
ters of galaxies as gravitational telescopes; A1689, A1835, and 
AC114. A1689 {z = 0.183) has the largest Einstein radius o f all 
massive lensing clusters, 9e ~ 50". lBroadhurst et al.l ( |2005[) and 
iLimousin et al l (l2007h performed a strong lensing analysis us- 
ing HST data and identified 115 images of 34 multiply lensed 
background galaxies in the redshift range 1 < z < 5.5 based on 
spectroscopic and photometric redshift estimates. A cluster mass 
model of AC114 (z = 0.312, Oe ~ 30") and several strongly 
lensed sources, i ncluding a 5 -image confi guration at z = 3.347 
were presented in lCampusano et al] (1200 lb . The mass model for 
A1835 (z = 0.253) yields an Einstein radius of 9e ~ 40" at 
high-z jRichard et al.Ll2006l) . 



4.1. Lensing magnification: tunnel vision 

To calculate the lensing magnification of the SN lightcurves, 
the public LENSTOOlfl software package was used. The code 
is specifically developed for modeling the mass distribution of 
galaxies and cluste rs in the strong and weak lensing regime 
jKneib et aLl Il996l). It u ses a Monte Carlo Markov Chain tech- 
nique (lJullo et al. n2007h to constrain the parameters of the clus- 
ter model using observational data of the background galaxies 
as input. The output can then be used to compute, e.g., the mag- 
nification and time delay function at any given position behind 
the cluster For A 1689 the mass model by lLimousin et al ] (120071) 
was used. Th e clust ers A1835 and AC114 were modeled as in 
iRichard et all (|2006|) . Figure |4] shows the average lensing mag- 
nification as a function of source redshift in the FOV of the 
ISAAC camera, which is 2.5 x 2.5 arcmin^ for the three clus- 
ter fields considered. We note that A 1689 seems to be the most 
promising gravitational lens for reaching the highest redshifts. 
For the 2007 observations of A 1689 centered on the cluster it- 
self, the magnification is on average, >2.5 mag for z^l in the 
ISAAC field of view. The 2003/2004 archival observations that 
we used were off'set from the cluster core, and therefore the av- 
erage magnification for these observations is lower, ~1.5 mag 
for z^l. For A1835, the average magnification is ~1 mag for 
Z>1. The width of the A1689 2003/2004 and A1835 curves indi- 
cate the slightly different pointings and effective FOV of these 



observations. AC114 has an average lensing magnification of 
~0.8 mag for z>l. 
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Fig. 4. Average lensing magnification versus redshift for the 
three observed clusters computed with LENSTOOL (see text) for 
the ISAAC/VLT field of view. Two curves are shown for A 1689, 
since the pointings in 2003/2004 (dashed red) and 2007 were 
different. The width of the curve for the A 1689 pointings (solid 
red) in 2003/2004 and the A1835 pointings (orange) indicate 
that these pointings do not cover the same area, but are slightly 
shifted between observations and the effective area observed is 
different (details can be found in Paper I). 



4.2. Monitoring time for SN surveys witli gravitational 
telescopes 

Because of flux conservation, large lensing magnifications result 
in small observed solid angle co. 

Therefore, unlike other SN searches, the effective solid angle 
CO in Eq. (|6]l is not constant with redshift when cluster fields are 
targeted. The light beam at any given redshift behind the cluster 
is magnified by a factor 



(10) 



(11) 



^ _ ]^Q-0.4 Ami„s 

at the expense of a smaller solid angle to of viewing 

where dto and 6loq represent infinitesimal solid angle elements, 
with and without lensing magnification. Therefore, the effective 
volume of a lensed SN search dVc can be measured by substi- 
tuting the expression into Eq. ^ 



dVr 



cdliz) _ ^ 

rW dz. 



H(z)(i+zy 



(12) 



www. oamp . fr/cosmology/lenstool 



Throughout this paper, A/wiens < for ^ > 1 . 
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The corresponding reduction in the source area as a function of 
redshift for the strongest (A1689) and weakest (AC114) lens in 
our survey are shown in Fig.|5] Although A1689 is the the most 




[deg] 



Fig. 5. Source plane area shrinkage behind our strongest gravita- 
tional telescope A1689 (top figure) for various redshifts; z =0.25 
(orange), 0.5 (yellow), 2.0 (magenta) and the weakest in the cur- 
rent survey ACl 14 (bottom figure) for various redshifts (z = 0.5 
(yellow), 2.0 (magenta), 3.0 (cyan). The utmost line (red) shows 
the effective FOV of the observations. 

promising gravitational lens, the figures illustrate that the rapid 
shrinkage in the solid angle with increasing redshift is also very 
strong. 

Thus, as shown in iGunnarsson & Goobad (l2003h . a gravi- 
tational lens does not always enhance the number of SN dis- 
coveries. However, it does increase the limiting redshift of a 
magnitude-limited survey. We exploit this effect to search for 
SNe at redshifts beyond those explored by "traditional" SN 
searches. 

Weaker lenses, such as ACl 14 and A1835, may not go as 
deep in redshift, but one may still find a comparable number of 
SNe as behind A 1689 (or even more), although these SNe would 
be found at somewhat lower average redshifts. 



' Gravity gives, gravity takes ! 



For the unlensed case, the monitoring time above threshold 
for a SN of type j, T is a function of the SN lightcurve, the de- 
tection efficiency, e, the extinction by dust, Amext, and the intrin- 
sic brightness Mj with the probability distribution P{Mj). With 
Atj being the lightcurve time period when the supernova is above 
the detection threshold, 

Tjiz, Amext) = e ■ J Af y(z, Mj + Am,,,) P{Mj)dMj. ( 1 3) 

We assume P{M) to be Gaussian, and to have the mean values 
and standard deviations listed in Table [T] 

Taking into account the lensing effect of the clusters, the 
monitoring time becomes 

r/z, Amext, Amiens) = 

€ ■ J Atjiz, Mj + Amext + Amiens) P{Mj)dMj, (14) 

keeping in mind that Amext > corresponds to dimmed SNe. 
Usually, Amiens < 0, and SNe will be magnified (although there 
are also areas of the field where the gravitational lens demagni- 
fies). 

Thus, the expected number of SNe for a given type j, using 
volumetric rates, is then given by 

Nj = Tj(z, Amext, Amiens) ■ ■ dVc, (15) 

J (1+z) 

where we assume a n overall Milky-Way-like dust extinction 
jCardelhet all 1 19891) with E{B - V) = 0.15 and Ry = 3.1. In 
this study, we have assumed a constant optical depth for all su- 
pernovae. This choice matches the assumptions used to derive 
the SFR estimates that we have used. 

5. Properties of background galaxies 

For each one of the three considered clusters, galaxy catalogs 
were compiled using archival optical and near-lR photometry. 
The different instruments and filters are listed in Table |2] 

The observed magnitudes in at least three bands, optical or 
near-IR, were used to derive a photometric r e dshift using the 
templa te-fitting technique (e.g., iGwvnl (1 1995b : iMobasher et al.l 
(ll996h V The software used is the c ode de veloped by the GOODS 
team as described in Dahlen et al. ( 120051 : 2009, in prep.) 

Figure |6] shows the distribution of galaxies in bins of Az = 
0.5 for the three clusters. Furthermore, the restframe UV-flux 
at 2800A of each resolved galaxy was computed and used as a 
tracer of the SFR. Similarly, the integrated SFR was calculated in 
each redshift bin by summing the resolved galaxies behind each 
cluster. Thus, we compute the expected number of SNe with two 
methods: 1) the volumetric rate taken from the literature, and 2) 
the measured star-formation rate of the resolved galaxies in the 
FOV. 

6. The ISAAC pilot survey 

During the spring of 2007, three clusters, A1689, ACl 14 and 
A1835, were monitored with the ISAAC instrument at VLT of 
approximately one-month intervals, as shown in Table [3] In to- 
tal, the data set consists of 20 hours of VLT time on target: 4.5 
hr, 5.85 hr and 10 hr for A1689, ACl 14 and A1835 respectively. 
The data is complemented by archival data (also listed in the ta- 
ble): 8.4 hr, 5.7 hr and 6.1 hr for A1689, ACl 14 and A1835. 
The survey filters in our 2007 program were chosen to match the 
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Filter 




Instrument/Camera 


Abell 1689 


F475W, F625W, F775W, F850LP 


HST/ACS 


FUOW, F160W 




HST/NICMOS 


Abell 1835 


V, R, I, 




CFHT/CFH12K 


F702W 




HST/WFPC2 


Z 




VLT/FORS2 


SZ, J, H, Ks 




VLT/ISAAC 


AC 114 


U 




CTIO 


B 




AAT/CCDl 


V 




ESO-NTT 


F702W, F814W 




HST/WFPC2 


J, H, Ks 




VLT/ISAAC 



450 
400 
350 
.2 300 
250 
^ 200 
^ 150 
100 
50 






A1689 ■ 




AC114' 




A1835 ■ 


! p.^.™ 









Fig. 6. Resolved galaxies vs redshift behind A1689, A1835 and 
AC 114. 



deepest reference images. Thus, the SN search was done using 
the 5Z-filter for A1689 and A1835 and /-band for AC114. A 
full description of the data reduction, SN search efficiency, and 
limiting magnitude is reported in Paper I. An average discov- 
ery depth at 90 % CL of 5 Z, 7 < 24 mag (Vega) was derived by 
Monte Carlo simulations in which artificial stars were added to 
the images. 

6.1. Expected event rate in the survey 

One of the most important aspects of the pilot survey is to ex- 
plore whether the use of gravitational telescopes significantly 
enhances the survey depth given the observational magnitude 
limit. In the upper panels of Figs. |7] [8] and |9] we explore the 
differential number of SNe expected for each one of the three 
clusters with lensing magnification. The lower panels of the fig- 
ures show the gain/loss due to the lensing compared to the same 
survey without lensing as a function of redshift. In particular, 
the lower panels indicate the redshift regions where the use of 
gravitational telescopes enhances the detection probability. 

As expected, the boost is most important for the fainter core- 
collapse supernovae. Type Ib/c and IIP in particular, where the 
detection efficiency is increased for z > 0.5. For the brighter 
SNe, such as Type la, it is only for z > 1 that a net gain is ex- 
pected. Thus, the foreground massive cluster, besides increasing 
the flux levels, serves as a high-z filter. 

Both AC114 and A1835 (to a somewhat lesser extent) pro- 
vide comparable total SN rates to A 1689 although A 1689 is 



Date 


Exposure 


Seeing 


90% detection 


Area^ 




[min] 


[arcsec] 


efficiency [mag] 


[arcmin^ 




Abell 1689 


-VLT/ISAAC 5 Z-band 




2003 02 09* 


159 


0.52 


24.28, transient 




2003 04 27 


43 


0.43 


transient 


3.70 


2004 01 13 


43 


0.52 


23.58, non-detect 


2004 02 14 


43 


0.58 


23.64, non-detect 




2003 01 16 


43 


0.58 


23.48 




2003 02 15 


43 


0.50 


23.60 


3.72 


2003 04 27 


86 


0.44 




2004 01 12 


43 


0.55 


23.64 




2007 04 08 


117 


0.64 


23.95 




2007 05 14/15 


117 


0.65 


23.95 


4.44 


2007 06 06 


39 


0.70 


23.15 






AC 114- 


- VLT/ISAAC /-band 




2002 08 20 


108 


0.49 


23.87 




2007 07 13" 


234 


0.43 


24.04 




2007 08 09 


117 


0.73 


23.79 


5.06 


2007 09 02 


117 


0.55 


23.83 




2007 09 28 


117 


0.46 


24.04 






Abell 1835 


-VLT/ISAAC 5 Z-band 




area 1 










2004 04 20 


231 


0.49 


24.06 


3.75 


2004 05 15 


135 


0.62 


24.06 


3.75 


2007 04 18 


117 


0.79 


23.80 


2.50 


2007 05 18 


78 


0.74 


23.83 


3.75 


2007 07 18 


117 


0.62 


23.80 


2.50 


area 2 










2007 04 18 


117 


0.79 


23.70 


1.57 


2007 05 14/18 


60 


0.80 


23.45 


2.13 


2007 07 18 


117 


0.62 


23.70 


1.57 



" - average of observations obtained on July 1 1,12,13 and 15. 

- average of observations obtained on February 5,1 1 and 15. 

- overlap region with other images to which the detection limit in 
column 4 applies. 



a much stronger lens than the other two clusters. This is be- 
cause, as already mentioned, the magnification is associated with 
a shrinkage in the effective volume element. 

For simplicity, we restricted this comparison in Figs. |7] to 
|9]to the volumetric rates estimates in lMannucci et al. 1 (120071) . In 
Fig.Hni the various model predictions for the three clusters com- 
bined are shown for each SN type separately. 

The expected total number of SNe in our survey is shown 
in Fig. [T T| The estimated rate s assume extinction by Galactic- 
like dust ICardenietalllI989h with an average color excess of 
E(B - V) - 0.15 and a total-to-selective extinction coefficient 
Ry - 3.1,i.e., Av - 0.46 mag. Since the lensing magnification is 
typically Amiens-5- 0.8, the impact from dust extinction accounts 
for less than a factor two decrease in the expected number of SN 
discoveries, compared to the results where dimming by dust is 
completely neglected. 

6.2. A transient candidate 

Paper I described the image subtractions used to search for tran- 
sient objects in our data set. Transients were sought in both the 
new images, using archival data as a reference, and for transients 
in the archival images using the period 79 (April-July 2007) data 
as a reference. The images were geometrically aligned and the 
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Fig. 7. Upper: Redshift distribution for the nu mber of SNe 
(for ea ch type) assuming the rates estimates by M annucci et al.l 
(l2007h for A 1689 in 5Z-band. Lower: Gain/loss of using A 1689 
as a lens compared to an equivalent survey without the lens for 
different redshifts. The crossing of the curves through the zero 
line indicates the redshift for which a transition to a net gain 
in SN discoveries is obtained due to the gravitational telescope. 
An average Milky-Way like extinction with E(B-V)=0.15 was 
assumed for both plots. 



Fig. 8. Upper: Redshift distribution of the expected num- 
ber of SNe (for each type) assuming the rates estimates by 
iMannucci et all (l2007h for A1835 in 5Z-band. Lower: Gain/loss 
of using A 1835 as a lens compared to an equivalent survey with- 
out the lens for different redshifts. The crossing of the curves 
through the zero line indicates the redshift for which a transi- 
tion to a net gain in SN discoveries is obtained due to the grav- 
itational telescope. An average Milky- Way-like extinction with 
E(B-V)=0.15 was assumed. 



Table 4. Transient candidate photometry from Paper I. 



Date 


Filter 


Magnitude (mag) 


2003-02-06 


/ 


24.09 ± 0.20 


2003-02-09/10 


z 


23.93 ± 0.08 


2003-02-26/27 


z 


23.94 ± 0.09 


2003-02-09 


SZ 


23.24 ± 0.08 


2003-04-12 


J 


23.61 ±0.15 


2003-04-27 


SZ 


23.73 ±0.16 



point-spread functions and the flux levels of the two images were 
matched prior to the pixel subtraction. 

In this process, one transient candidate was found in the 
A 1689 archival images in the S'Z-band. One /-band data point of 
A 1689 and two z-band data points were observed with FORS2 
and another in 7-band with ISAAC at VLT during the time that 
the transient was bright. We used HAWK-I /-band images from 
our program in period 81 (July 2008) as a reference to obtain a 
measurement of the transient flux in that band. Additional ref- 
erence data from F0RS2 and HST/ACS were available for the 
optical bands. These were used to measure the flux in the re- 
gion of the transient after it had faded. The transient photometry 
is summarized in Table |4] The location of the transient and the 
lensing magnification map is shown in Fig. [12] 

All available SN templates and a grid of redshifts (z - [0, 3]) 
and reddening parameters were tested (allowing for an intrin- 
sic variation in the brightness) and the best fit was found for 
a Type IIP te mplate based on lightcurves of SN2001cy from 
iPoznanski et al. (2009), redshifted to zs.n - 0.59. Moreover, 
the best fit of the transient colors was found by assuming that 
the SN is highly reddened, with a low total-to-selective ex- 



tinction ratio (Ay=1.27, Rv=l-5), as shown in Fig. \T3\ We 
note that low values of Ry, although not seen in the Milky 
Way, were reported for extinction of quasars dWang et all 120041: 
lOstman et all |2008[) and shown to be v ery common along SN 
lines of sight (iNobili & Goob^ l2008l) . poss ibly as a result 
of multiple scattering by circumstellar dust jGoobait l2008h . 
Another possibility is that the intrinsic colors of the SN candi- 
date differ significantly from SN2001cy, in which case a bias 
could be introduced in the K-corrections. H owever, a recent 
study of 40 low-z Type IIP SN lightcurves dPoznanski et all 
l2009h found a low average value of the total to selective extinc- 
tion ratio, - 1.5 +0.5, in excellent agreement with the best-fit 
solution for our SN candidate. 

For the nearest galaxy, at 0.5" projected distance, a photo- 
metric redshift Zgai = 0.60 + 0.15 is derived, as shown in Fig. [14] 
The closest galaxy, 1 .2" away, has a photometric redshift that 
also peaks at z = 0.6. The distance to the galaxy cores is thus 3.3 
and 7.7 kpc, respectively. 

It should be noted that at z = 0.6 because of the dust ex- 
tinction the transient (/ ~ 25 mag, unlensed) would not have 
been detected in our survey without the magnification power of 
the cluster, 1.4 mag. Taking into account the lensing magnifi- 
cation and the assumed host galaxy extinction, we find that the 
best fit shown in Fig. [T3] corresponds to an absolute magnitude 
Mv - -17.6 + 0.3, in good agreement with the assumptions in 
Table [1] It is striking that the tentative identification of the tran- 
sient redshift and type match well the expectations for the survey 
in terms of SN subtype and redshift, as shown in Fig. [7] 

We now address the alternative possibility that the tran- 
sient is at the cluster redshift. A Type la supernova more than 
70 days past lightcurve maximum co uld potentially matc h the 
observed brightness of the transient. ISharon et al.l (l2007h esti- 
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AC 114 




Fig. 9. Upper: Redshift distribution of the expected num- 
ber of SNe (for each type) assuming the rates estimates by 
iMannucci et alJ (l2007h for AC114 in 7-band. Lower: Gain/loss 
of using ACl 14 as a lens compared to an equivalent survey with- 
out the lens for different redshifts. The crossing of the curves 
through the zero line indicates the redshift for which a transi- 
tion to a net gain in SN discoveries is obtained due to the grav- 
itational telescope. An average Milky-Way-like extinction with 
E(B-V)=0.15 was assumed. 



mated the rate of Type la supernovae in intermediate-z clusters 
to ~ 0.3 SNu. For the integrated galaxy luminosity in A 1689 of 
Lb ~ 1.25 • IO'^Lq and a monitoring time of ~ 100 days, about 

0. 1 Type la supernovae were expected in A1689 in our dataset, 

1. e., a non-negligible possibility. However, as discussed in paper 
1, the decline rate of SNla lighturves at late times is 0.01-0.02 
mag/day, which is inconsistent with the transient lightcurve. 

To conclude, we find that the match to the photomet- 
ric redshifts of the potential host galaxy and the supernova 
lightcurves, along with the fitted peak magnitude being con- 
sistent with a Type IIP supernova (reddened by dust with Ry 
similar to that found fo r the nearby sample of Type IIP SNe in 
dPoznanski et al.L l2009l) ) to represent the most compelling fit to 
the transient behind A 1689. 



7. Implications for future near-IR surveys 

The pilot survey was completed with the ISAAC instrument 
at VLT, which has a FOV of 2.5'x2.5' and a threshold of 
~24 mag (Vega) for SZ and /-bands and relatively few observa- 
tions. We briefly discuss the feasibility of building up lightcurves 
of lensed SNe behind clusters of galaxies for surveys with 8- 
meter class ground-based telescopes and large FOV near-IR in- 
struments, such as HAWK-1 at VLT or MOIRCS at Subaru. We 
consider a five-year "rolling" search survey, with imaging spaced 
at intrevals of 30 days. 

The lensing model of A 1689 is used as our baseline for esti- 
mating the number of SNe behind a massive cluster as a function 
of redshift. Thus, the results below apply to the most massive 
clusters, M ~ lO'^M©. We also consider a survey period of five 
years since this is optimal for the discovery of multiple images 
(Sect. 17. It . In practice, several clusters would have to be ob- 
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Fig. 10. Redshift distribution (including an overall reddening of 
E(B - y) = 0.15 and Ry = 3.1) for the number of SNe for the 
various model predictions for the three clusters combined for 
each SN type separately. 



served to correspond to 'five A1689 years' since that particular 
field is behind the Sun part of the year 

During period 81 (P81, July 2008), our team started a survey 
targeting lensed SNe behind A 1689 using the increased sensi- 
tivity and FOV (* 7.5'x7.5') of HAWK-1 on VLT. The Umited 
availability of the instrument in P8 1 , only a few, closely spaced 
observations were completed. Although the observations were 
not suited to transient searches, the dataset could be used to es- 
timate the depth of the point-source search in /-band that was 
continued during P82 (early 2009) to ~ 24.65 mag (Vega) for 
90% detection efliciency, i.e., about 0.65 mag deeper than the 
survey done with ISAAC. 

We now examine the feasibility of SN detection with 
HAWK-1. We use LENSTOOL to calculate the lensing magnifi- 
cation map of A 1689 for the larger FOV. Although the magni- 
fication decreases with distance from the cluster core, reaching 
Amiens ~ -0.25 at the edges of the FOV, the impact of lensing 
remains very important for detecting distant supernovae. In the 
upper panel of Fig. [15] the differential number of SNe is shown 
for the various types of SNe. The lower part of the figure shows 
the gain and loss of lensing compared to the same survey with- 
out lensing as a function of redshift. As expected, fewer SNe are 
found at low redshifts, while the survey depth is significantly in- 



10 



A. Goobar et ah: Near-IR search for lensed supernovae behind galaxy clusters 



A1689 + A1835 + AC114 



la ■ IIP □ lln 
HN □ Ibc ■ IIL 





M07 



KN08 



D04/D08 



SFR 



A1689 + A1835 + AC114 




N06 



SB05 



TAUl 



Fig. 11. Number of expected SNe for the observations given in 
Table |3] for (1) the SN rates shown in Figs. |2] and |3] and (2) 
the SN rates derived from SFR (denoted by SFR) in resolved 
galaxies described in Sect. 13.21 
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Fig. 12. LENSTOOL lensing magnification map of A 1689 
based on mass model in [Limousin et al.l (l2007h superim- 
posed on top of a HAWK-I 7-band image from our P81 
program (July 2008). The lensing contours for a source at 
z - 0.6 for Amiens=(-0.5, -1, -2, -5) mag shown in (or- 
ange,red,yellow,green). The arrow points to the position of the 
transient, where a magnification of 1.4+0.3 mag is expected. 



creased. The integrated number of SNe for the various models is 
shown in Fig. [16] 

For HAWK-1 (and the A 1689 mass model), we expect on 
the order of 40-70 SNe (depending on the underlying rate esti- 
mate for the various SN types) of which about a dozen will be at 
z > 1-5. In Fig.[T71 the potential of the suggested rolling search 
for generating lensed supernova lightcurves in the redshift range 
(1 < z < 2) is shown. The repeated images are used both to dis- 
cover new supernovae and to build up lightcurves for earlier dis- 
coveries. The supernova types and redshift distribution matches 
the differential rates in Fig. [15] 

Surveys for lensed supernovae with space instruments would 
complementthe ground-based approaches since even higher red- 
shifts could be reached due to the deeper point source sensitivity. 

7. 7. Multiple SN images 

When looking through a gravitational lens, multiple images of 
the same source image can be observed. This is also true for 
SNe that, due to strong lensing, can potentially be detected to 
very high redshifts. About one in a hundred SNe behind A1689 
in the HAWK-1 FOV would have multiple images with time sep- 
arations of between weeks and a few years. Thus, about 0.5-1.0 
SNe with multiple images are expected in a 5 year survey with 
HAWK-1. 

Figure[T8]indicates the fraction of the source areas with mul- 
tiply lensed SNe that can be observed as a function of survey 
time for two different redshifts. For z - 1, all SN types show the 
same behavior and given a sufficiently large survey time, at least 
two (or more) images of the SN could be observed, regardless 
of its type. At higher redshifts (z - 2), given a sufficiently large 



survey time, most of the brighter SNe (la and lln) and about half 
of the llLbright and HN will be observable and have at least two 
(or more) images. The other SN types (IIL, Ib/c, and IIP) will 
be too faint - even with magnification - to be observed. For all 
considered redshifts, a 5 year survey (or longer) is optimal for 
discovering multiple images of SNe behind clusters. 

Detecting these rare events could provide important con- 
straints on the Hubble constant with the time-delay technique 
dRefsdall 1 19641) as well testing dark mat ter and energy prop- 
erties in an unexplored re dshift range (iGoobaret al.L 120021: 
M ortsell & Sunessonll2006l) . Models of lens systems are in gen- 
eral uncertain due to the possibility to rescale the density dis- 
tribution of the lens and add a circularly symmetric density 
mass-sheet, while preserving the observed i mage configuration 
dGorenstein et al.L 119881: iLiesenborgs et all |2008|) . This mass- 
sheet degeneracy can be broken if the absolute magnification 
of the lens is known. Since SNla have a very tight dispersion 
in brightness, these lens systems would constitute an ideal sam- 
ple for minimizing three major systematic uncertainties in the 
estimates of Hq using the time-delay technique: accurate time 
estimates from supernova lightcurves, elimination of the mass- 
sheet degeneracy, and accurate lens models because of the large 
number of lensed background galaxies, as discussed in Paper 111. 

8. Conclusions 

Powerful gravitational telescopes in the form of massive galaxy 
clusters provide unique opportunities to discover transient ob- 
jects such as SNe at redshifts beyond what can be reached with 
current telescopes. The lensing magnification fi corresponds to 
a gain factor in exposure length, fr, while at the same time the 
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Fig. 17. Simulated lightcurve sampling for supernovae in the redshift range 1 < z < 2 in a 5-year monthly survey of A1689 and 
AC114 with HAWK-1. The top panel shows the expected lightcurves sampling SNla. Core-collapse SNe are shown in the bottom 
panel. The vertical axis shows the unlensed magnitudes. The symbol colors indicate the SN redshift: blue (z - [1.0, 1.2]), cyan 
(z= [1.2, 1.4]), green (z = [1.4, 1.6]), orange (z = [1.6, 1.8]), and red (z = [1.8,2.0]). 



solid angle at the source planes shrinks by a factor ju for source 
redshifts higher than the cluster redshift. The net gain/loss of 
searching for supernovae behind massive clusters is therefore a 
non-trivial combination of FOV, limiting depth, and supernova 
luminosity functions. In general, the lens works as a magnifying 
glass and high-z filter, i.e., by reducing the number of detections 
of bright/close supernovae and enhancing the detections of dis- 
tant/faint objects. Thus, to be successful, a SN survey behind 
galaxy clusters needs to be optimized. Clearly, for extremely 
sensitive (e.g JWST or ELT) or very large FOV instruments, the 
positive impact of the lensing cluster may be negligible, at least 
for the brightest types of supernovae. The net benefit of exploit- 
ing the suggested approach will ultimately depend on the rate 
and intrinsic brightness of the various types of SNe at redshifts 
beyond what is currently known. For Type la supernovae, an im- 
portant parameter determining the rates beyond z - 1 .5 is the 
delay time, t. By increasing the redshift sensitivity beyond that 
achieved by "standard" surveys, we may significantly improve 
our understanding of SNla progenitors. Similarly, little is known 
about the dimming of supernovae by dust at very high redshifts. 
The combination of a longer wavelength-survey and higher sen- 
sitivity to fainter high-z SNe could thus lead to detections of a 
different population of objects. 

A combined 40-hour dataset involving archival ISAAC data 
and new observations obtained in 2007 for three very massive 
clusters (A1689, A1835, and ACl 14) was used to determine the 
feasibility of discovering lensed core-collapse and Type la SNe. 
Considering the monitoring time available, the area surveyed, 
the lensing magnification, and the survey magnitud limit, rate es- 
timates of the various SN subtypes considered were calculated. 
Synthetic lightcurves of SNe and several models of the volumet- 
ric Type la and core-collapse SN rates as a function of redshift 



were used, all consistently predicting a Poisson mean value for 
the expected number of SNe in the survey of between Nsn=0.8 
and 1.6 for all SNe. One transient object was found behind 
A1689 on a galaxy with photometric redshift Zgai - 0.6 + 0.15, 
the most probable redshift for SN detection in the ISAAC/VLT 
survey. The lightcurve is consistent with being a reddened Type 
IIP supernova at zsn = 0.59 + 0.05. At the position and redshift 
of the transient, the lensing model predicts 1.4 magnitudes of 
magnification. 

Because of the recent deployment of large and sensitive 
near-lR cameras, such as HAWK-1 at VLT, the search for the 
highest redshift SNe can now be moved to longer wavelengths, 
thus avoiding the difficulties involved with restframe UV ob- 
servations, and extending the potential for supernova discov- 
eries, especially Type la supernovae, beyond z > 2. A feasi- 
bility study of the potential to build up lightcurves of lensed 
SNe with larger and deeper surveys shows that this is a very 
exciting path for new discoveries. The equivalent of a five-year 
monthly survey of a single very massive cluster with the HAWK- 
1 camera at VLT would yield 40 - 70 lensed SNe, most of them 
with good lightcurve sampling. Thus, a dedicated multi-year 
NIR rolling search targeting several massive clusters would lead 
to a high rate of very high-z SN discoveries, thus making this 
appro ach complementar y to deep optical space-based SN sur- 
veys ( iRiess et al.L |2007|) as w ell large field-of-view optical SN 
searches, e.g., dPoznanski et al.. 2007.). 



Although very rare, multiple images of strongly lensed SNe 
are within reach of such a survey and could offer potentially ex- 
citing tests of cosmological parameters as well as improvements 
to the cluster mass modeling. 
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Fig. 13. Transient candidate photometry (also listed in Table |4]l 
plotted on top of redshifted (zsn = 0.59 ) lightcurves of the ver y 
well observed Type IIP SN SN2001cv |Poznanski et al] (|2009|) . 
The nearby SN lightcurves were K-corrected using IIP spectral 



templates derived from spectra and photom etry of SN2001c\ 

ricSOHi 
1.5 yielding 



Ley 

and reddened following the extinction law in iFitzpatrickl ( 1 1 9991) 
with the parameters E{B - V) = 0.85 and Ry - 
= 3.5 for 6 data points and 4 free parameters. 



Fig. 15. Upper: Redshift distribution of SN discoveries in a 
5-year survey behind a very massive cluster (model of A 1689 
used). Lower: Gain/loss of using a A1689-like cluster as a lens 
compared to an equivalent survey without the lens for different 
redshifts. The crossing of the curves through the zero line in- 
dicates the redshift for which a transition to a net gain in SN 
discoveries is obtained due to the gravitational telescope. An av- 
erage MiUcy- Way-like extinction with E{B -V) - 0.15 was as- 
sumed together with SN rates from M07. 




Fig. 14. Probability distribution function for the photometric red- 
shift (seven bands) of the galaxy closest (0.5 arcseconds) to 
the transient, the SN photo-z assuming a IIP template matching 
SN2001cy and the product of the two, yielding z - 0.59 + 0.05. 



Acknowledg merits 

We would like to thank Peter Nugent for providing lightcurve 
and spectral templates used in this analysis. Filippo Mannucci 
is also thanked for making his SN rate predictions available to 
us. We are also grateful to Dovi Poznanski for providing us with 
lightcurves and spectra of SN2001cy and to Avishay Gal- Yam 
for comments on an earlier draft. KP gratefully acknowledges 
support from the Wenner-Gren Foundation. AG, VS and SN 
acknowledge support from the Gustafsson foundation. V.S. ac- 



knowledges financial support from the Fundagao para a Ciencia 
e a Tecnologia. AG and EM acknowledge financial support from 
the Swedish Research Council. JPK thanks CNRS for support 
as well as the French-Israeli council for Research, Science and 
Technology Cooperation. 

References 

Astier, P. et al. 2006, A&A, 447, 31 
Broadhurst, T. et al. 2005, ApJ, 621, 53 

Campusano, L. E., Pello, R., Kneib, J. -P., Le Borgne, J.-F., Fort, B., Ellis, R., 

Mellier, Y, & Small, I. 2001, A&A, 378, 394 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 

Dahlen, T., Mobasher, B., Dickinson, M., Ferguson, H. C, Giavalisco, M., 

Kretchmer, C, & Ravindranath, S. 2007, ApJ, 654, 172 
Dahlen, T., Mobasher, B., Somerville, R. S., Moustakas, L. A., Dickinson, M., 

Ferguson, H. C, & Giavalisco, M. 2005, ApJ, 631, 126 
Dahlen, T., Strolger, L.-G., & Riess, A. G. 2008, ApJ, 681, 462 
Dahlen, T. et al. 2004, ApJ, 613, 189 
Filippenko, A. V. 1997, ARA&A, 35, 309 
Fitzpatrick, E. L. 1999, PASP, 111, 63 
Frieman, J. A. et al. 2008, AJ, 135, 338 

Gal- Yam, A., Maoz, D., & Sharon, K. 2002, MNRAS, 332, 37 
Giavalisco, M. et al. 2004, ApJ, 600, L103 

Gnedin, N. Y, Kravtsov, A. V., & Chen, H.-W. 2008, ApJ, 672, 765 
Goobar, A. 2008, ApJ, 686, L103 

Goobar, A., Mortsell, E., Amanullah, R., & Nugent, P 2002, A&A, 393, 25 
Gorenstein, M. V., Shapiro, 1. 1., & Falco, E. E. 1988, ApJ, 327, 693 
Gunnarsson, C, & Goobar, A. 2003, A&A, 405, 859 

Gwyn, S. D. J. 1995, Master's thesis, MS Thesis, University of Victoria (1995) 
Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 

Julio, E., Kneib, J. -P., Limousin, M., Eh'asdottir, A., Marshall, P. J., & Verdugo, 

T. 2007, New Journal of Physics, 9, 447 
Kim, A., Goobar, A., & Perlmutter, S. 1996, PASR 108, 190 



A. Goobar et al. : Near-IR search for lensed supernovae behind galaxy clusters 



13 



HAWK-I - 5 years - A 1689 



la □ IIP □ Iln 
HN □ Ibc ■ IIL 



E(B-V) = 0.15 



M07 



KN08 



D04/D08 



HAWK-I - 5 years - A 1689 




N06 



SB05 



TAUl 



HAWK-I - 5 years - A 1689 - z > 1.5 



■ la □ IIP □ Iln 
□ HN □ Ibc ■ IIL 



M07 



KN08 



D04/D08 



HAWK-I - 5 years - A 1689 - z > 1.5 




N06 



SB05 



TAUl 



Fig. 16. {Top) Number of SNe expected in a 5 year monthly 
survey of one very massive, A1689-like cluster, with HAWK-I. 
(Bottom) Number of SNe with z> 1.5. 
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Fig. 18. Fraction of SNe with multiple images that are observ- 
able as a function of survey duration. 



Kneib, J.-P., Ellis, R. S., Santos, M. R., & Richai'd, J. 2004, ApJ, 607, 697 
Kneib, J.-R, Ellis, R. S., Small, 1., Couch, W. J., & Sharpies, R. M. 1996, ApJ, 
471, 643 

Kobayashi, C, & Nomoto, K. 2008, ArXiv:0801.0215 

Kobayashi, C, Tsujimoto, T., Nomoto, K., Hachisu, I., & Kato, M. 1998, ApJ, 
503, L155+ 

Kolatt, T. S., & Bartelmann, M. 1998, MNRAS, 296, 763 
Kovner, I., & Paczynski, B. 1988, ApJ, 335, L9 
Kuznetsova, N. et al. 2008, ApJ, 673, 981 

Leibundgut, B. 2008, General Relativity and Gravitation, 40, 221 
Lien, A., & Fields, B. D. 2009, Journal of Cosmology and Astro-Particle Physics, 
1,47 

Liesenborgs, J., de Rijcke, S., Dejonghe, H., & Bekaert, P 2008, MNRAS, 386, 
307 

Limousin, M. et al. 2007, ApJ, 668, 643 

Mannucci, E, Delia Valle, M., & Panagia, N. 2007, MNRAS, 377, 1229 
Mannucci, E, Delia Valle, M., Panagia, N., Cappellaro, E., Cresci, G., Maiolino, 

R., Petrosian, A., & Turatto, M. 2005, A&A, 433, 807 
Miknaitis, G. et al. 2007, ApJ, 666, 674 

Mobasher, B., Rowan-Robinson, M., Georgakakis, A., & Eaton, N. 1996, 
MNRAS, 282, L7 

Mortsell, E., & Sunesson, C. 2006, Journal of Cosmology and Astro-Pai'ticle 

Physics, 1, 12 
Neill, J. D. et al. 2006, AJ, 132, 1126 
Nobih, S., & Goobai-, A. 2008, A&A, 487, 19 
Ostman, L., Goobai', A., & Mortsell, E. 2008, A&A, 485, 403 
Poznanski, D. et al. 2009, ApJ, 694, 1067 
. 2007, MNRAS, 382, 1 169 

Pritchet, C. J., Howell, D. A., & Sullivan, M. 2008, ApJ, 683, L25 
Refsdal, S. 1964, MNRAS, 128, 307 

Richard, J., Pello, R., Schaerer, D., Le Borgne, J.-E, & Kneib, J.-P 2006, A&A, 
456, 861 

Richardson, D., Branch, D., Casebeer, D., Millard, J., Thomas, R. C, & Baron, 

E. 2002, AJ, 123,745 
Riess, A. G. et al. 2007, ApJ, 659, 98 
Salpeter, E. E. 1955, ApJ, 121, 161 
Scannapieco, E., & Bildsten, L. 2005, ApJ, 629, L85 

Seitz, S., Saglia, R. P, Bender, R., Hopp, U., Belloni, P, & Ziegler, B. 1998, 
MNRAS, 298, 945 

Sharon, K., Gal- Yam, A., Maoz, D., Filippenko, A. V, & Guhathakurta, P 2007, 
ApJ, 660, 1165 

Smartt, S. J., Eldridge, J. J., Crockett, R. M., & Maund, J. R. 2009, MNRAS, 
395, 1409 

Stanishev, V. et al. 2009, ArXiv e-prints, 0908.4176 
Strolger, L.-G. et al. 2004, ApJ, 613, 200 



14 



A. Goobar et ah: Near-IR search for lensed supernovae behind galaxy clusters 



Sullivan, M., Ellis, R., Nugent, P., Small, I., & Madau, P. 2000, MNRAS, 319, 
549 

Sullivan, M. et al. 2006, ApJ, 648, 868 

Wang, J., Hall, R B., Ge, J., Li, A., & Schneider, D. P 2004, ApJ, 609, 589 



